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Abstract 

Wc study the motion of phase interfaces in a diffusive lattice equation with bistable nonlinearity 
and derive a free boundary problem with hysteresis to describe the macroscopic evolution in the 
parabolic scaling limit. 

The first part of the paper deals with general bistable nonlincarities and is restricted to nu- 
merical experiments and heuristic arguments. We discuss the formation of macroscopic data and 
present numerical evidence for pinning, depinning, and annihilation of interfaces. Afterwards 
we identify a generalized Stefan condition along with a hysteretic flow rule that characterize the 
dynamics of both standing and moving interfaces. 

In the second part, we rigorously justify the limit dynamics for single-interface data and a 

special piecewise affine nonlinearity. We prove persistence of such data, derive upper bounds for 

the macroscopic interface speed, and show that the macroscopic limit can indeed be described by 

the free boundary problem. The fundamental ingredient to our proofs is a representation formula 

0^ that links the solutions of the nonlinear lattice to the discrete heat kernel and enables us to derive 

.^^ macroscopic compactness results in the space of continuous functions. 
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1 Introduction 

Discrete forward-backward diffusion equations appear in many different applications such as edge- 
detection in digital images [PM90], models for population dynamics based upon random walks on 
lattices [HPO04], or phase transition problems with supercooling and superheating [E1185]. In all 
applications it is a common major problem to understand how the backward-parabolic regions affect 
the dynamics on large scales. In this paper, we study the diffusion lattice 



for j £ Z, t > 0, or equivalently 



Uj(t) = A$'(uj(t)) 



Wj(t) = V_$'(V+w/j(i)). 



XD 



X2) 



Here V_, V+ are the left- and right-sided discrete difference operators, A denotes the standard 
discrete Laplacian Apj = Pj+i — 2pj + Pj-i, and Uj, Wj are connected via Uj = V+uij. As illustrated 
in Figure 1.1, we always suppose that < I >/ the derivative of a double- well potential $, so it consists of 
two stable branches that enclose an unstable one. 




Figure 1.1: Left. Bistable derivative of a general double- well potential $. The increasing and decreasing 
branches of <F are called stable and unstable, respectively, while spinodal region refers to the interval [u*,u*] 
on which <£>' is decreasing. Right. Piecewise affine derivative of the degenerate double-well potential that is 
studied in §3 and corresponds to u# = —2, u* = u* = 0, vF = +2. 



Our goal is to characterize the effective dynamics of (1.1) in the parabolic scaling limit. Interpreting 
t > and j £ Z as the microscopic variables, we introduce the macroscopic time r > and space 
£GKby 

r = eh, £ = ej, (1.3) 

where e > is a small scaling parameter. We formally identify 



Ujit) = U(e 2 t, ej) 



(1.4) 



and aim to describe the evolution of U in the limit e 
the microscopic dynamics (1.1) into 



0. The scaling (1.3) and (1.4) transforms 



d T U = A £ &(U), 



(1.5) 



where A e is the standard finite difference approximation of dj on eZ, so the naive continuum limit 
as e — > reads 

d T U = dl<5>'{U). (1.6) 

This PDE, however, is ill-posed due to the unstable branches of $' and can therefore not determine 
the macroscopic limit of (1.1) completely. Actually, the lattice can be viewed as a regularization 
of (1.6) that accounts for small scale effects and provides in the limit e — >■ additional dynamical 
information such as laws for the motion of phase interfaces or even measure- valued solutions. Other 
notable regularizations of (1.6) are the Cahn-Hilliard equation 



d T u = dl$'(u) - e 2 dp, 



(1.7) 



which has enjoyed a lot of attention over the last decades, and the viscous approximation 

{l-e 2 dl)d T U = dl&(U), (1.8) 

studied in [NCP91, Plo94, EP04]. As discussed below, the available results indicate that the macro- 
scopic limits of the scaled lattice equation (1.5) and the viscous approximation (1.8) are identical 
but different from the limit of the Cahn-Hilliard model. In a formal way this can be understood by 
expanding the spatial operators in powers of e: Equations (1.1) and (1.8) share (up to a redefinition 
of e) the same leading order terms according to 

A £ P=(dl+ £ ^dl + 0{dl))p, (l-e%)- 1 dlP=(dl + e 2 dlP + 0(dl))P, 

with P being shorthand for $'(£/), while (1.7) replaces the right hand side in (1.6) by d 2 P — e 2 djU. 
For e > 0, however, the rescaled lattice and the viscous approximation are different and it remains 
open whether there exists a unified theory that is capable of describing the limit e — >■ for both 
models. 

A key feature of any regularization of (1.6) with double- well potential $ are phase interfaces, which 
evolve according to certain jump conditions and separate regions where U attains values in different 
stable regions (phases). Other types of nonconvex potentials give rise to different phenomena such 
as coarsening of localized spikes, see for instance [EG09]. 

Numerical simulations of (1.1) as performed in §2 with a generic double-well potential provide 
evidence for the existence of two different types of phase interfaces. Type-I interfaces correspond 
to piecewise smooth functions U and separate regions where U takes values in either one of the 
phases U < u* and U > u* , where [u*,u*] is the spinodal interval. Type-II interfaces, however, 
are related to measure-valued solutions of (1.6) and model a phase mixture on at least one side of 
the interface. For both types, a phase interface can have a fixed position or move, depending on 
the behavior of P = $'([/) near the interface. While P is smooth across a standing interface and 
takes values in [p*,p*] = [&(u*), $'(«*)], a moving interface is driven by a jump in d^P but requires 
either P = p* or P = p* subject to the propagation direction. In the macroscopic limit we therefore 
find hysteretic behavior in the sense that fronts moving into different phases comply with different 
constraints, see Figure 1.2. Further intriguing properties of the macroscopic lattice dynamics are 
sketched in Figure 1.3. Driven by the bulk diffusion, a standing interface can suddenly start to move 
(depinning) and a moving interface can eventually come to rest (pinning). Moreover, two interfaces 
can disappear after a collision (annihilation). 

Assuming that the lattice data pj = &(itj) converge as e — > to a sufficiently regular function P 
and that any phase interface is of type I, the dynamics in the parabolic scaling limit can described 
by combining bulk diffusion via (1.6) with the generalized Stefan condition 

|[P]|=0, ^\[U}\ + \[dsP]\=0 (1.9) 



dr 



and the hysteretic flow rule 



P=p, for ^-W]>0, P = p* for ^H[[/]<0, (1.10) 

dr dr 

where £*(r) is the position of an interface and [•] denotes the jump across this interface. These 
conditions have also been proposed in [EP04] to model the propagation of type-I interfaces in the 
limit of the viscous approximation and are naturally related to the notion of entropy solutions, see 
also [MTT09] and the discussion below. For the Cahn-Hilliard equation (1.7), the parabolic scaling 
limit does not imply any hysteresis. Here each interface corresponds to P = and evolves therefore 
according to the classical Stefan condition, see for instance [BBMN12] for a rigorous proof. We also 
note that there exists at least one other macroscopic limit for (1.7), which is, however, not related 
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Figure 1.2: Le/t. Cartoon of the macroscopic hysteresis for the potentials from Figure 1.1. The arrows 
indicate the temporal jump of U when it undergoes a phase transition at a fixed position £. In particular, 
P = p* holds at any interface that moves into the phase U < u* , whereas propagation into U > u* requires 
P =p*. The dashed lines represent standing interfaces, at which P takes values in [p*,p*]. 
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Figure 1.3: Cartoon of three macroscopic type-I interfaces. The first interface (moving) and the second one 
(standing) eventually collide with each other and disappear (annihilation). The third one is initially at rest, 
starts to move at a later time( depinning) and stops eventually again (pinning). 



to the parabolic scaling (1.3): In the regime of almost vanishing bulk diffusion, interfaces move and 
merge on a much slower time scale which is exponentially small in e [ABF91, BH92]. 

In the case of hysteretic interface motion, there seems to be no rigorous result - neither for the 
lattice nor the viscous approximation - that derives (1.9) and (1.10) rigorously from the dynamics for 
e > 0. Previous results for the lattices (1.1) or (1.2) are either restricted to standing interfaces, see 
[GNU] and [BGN13] for type-I and type-II interfaces, respectively, or do not capture the dynamics 
of moving interfaces completely, e. g. [BNP06]. 

As a first step towards a mathematical justification of the macroscopic evolution laws for type-I 
interfaces, we study in §3 the special case of 



$(u) 



2 mm 



{(u-l) 2 , (u + l) 2 }. 



(1.11) 



At the cost of being discontinuous at u = 0, the derivative $' of (1.11) has two advantages over a 
generic bistable function. First, the nonlinearity in (1.1) is piecewise afnne and second, the spinodal 
region has shrunk to the point u = 0, see the right panel of Figure 1.1. These properties simplify the 
dynamical system (1.1) significantly and enable us to represent solutions to the nonlinear lattice by 
a summation formula that involves delayed and shifted versions of the discrete heat kernel. 

Due to the degenerate nature of (1.11), it is not our intention to identify the most general 
class of admissible initial data for which the macroscopic limit can be described by a free boundary 
problem. On the contrary, in order to keep the presentation as simple as possible, we restrict our 
considerations in §3 to initial data that produce a single type-I interface which cannot change its 
direction of propagation. Our main findings are formulated in Theorems 3.2, 3.16, and 3.18, and can 
informally be summarized as follows. 

Main result. The lattice (1.1) with (1.11) has the following properties. 

1. Microscopic single-interface solutions: Type-I interfaces are naturally related to a class of lat- 
tice states that is invariant under the dynamics. Specifically, intervals of linear diffusion are 



interrupted by an increasing sequence of phase transition times (t* k ) k>k such that u k switches 
from negative to positive sign at t = t* k . 

2. Macroscopic evolution: When starting with macroscopic single-interface initial data, the limit 
e — > can be characterized as follows. 

(a) Convergence: The lattice data Uj converge in a strong sense to a function U that is 
smooth outside of an interface curve r i-> £*(t), where £* is Lipschitz continuous and 
nondecreasing. Moreover, the function P = U — sgn U is continuous across the interface. 

(b) Limit dynamics: For almost all times r > we have either 

P(r,C(r)) €[-!,+!] and ^t(r) = 
or 

P(r,C(r)) = +l and Af( r ) > , 

OlT 

and P solves the linear heat equation outside of the interface. 

(c) Uniqueness: U and £* are uniquely determined by the macroscopic initial data and the 
limit model. 

The rest of the paper is organized as follows. In §2.1 we employ the gradient flow structure of 
(1.1) to describe the formation of macroscopic data on a heuristic level. Afterwards we report on 
our numerical investigations for general double- well potentials. We present several examples for the 
macroscopic motion of type-I and type-II interfaces in §2.2, and proceed in §2.3 with discussing the 
key features of the microscopic dynamics near moving interfaces. These are: sequentiality of phase 
transitions, small-scale fluctuations, and existence of multiple time scales. Finally, in §2.4 we give 
a more detailed description of the macroscopic limit model for type-I interfaces and interpret the 
hysteretic flow rule in terms of entropy inequalities. 

§3 contains our analytical results for the special case (1.11). We first employ ODE arguments 
in §3.1 in order to prove the persistence of single-interface data. In particular, in Corollary 3.7 
we establish the aforementioned representation formula. In §3.2 we then introduce the concept of 
macroscopic single-interface initial data and derive upper bounds for the macroscopic interface speed 
from the properties of the discrete heat kernel, see also Appendix A. The main technical work is done 
in §3.3, were we establish macroscopic compactness results for £* and P in the spaces of Lipschitz 
and Holder continuous functions, respectively. In §3.4 we finally pass to the limit e — > 0. To this 
end, we first justify the limit model along subsequences, and obtain afterwards both uniqueness and 
convergence by adapting some arguments from the theory of free boundary problems with hysteresis 
operators. 

2 Heuristic arguments and numerical simulations 

In this section we employ heuristic arguments as well as numerical simulations in order to gain a 
qualitative understanding of the key dynamical features of the nonlinear lattice diffusion with double- 
well potential. In particular, we discuss (i) the underlying gradient flow structure and the formation 
of macroscopic data during a fast initial transient regime, (ii) the microscopic and macroscopic 
dynamics of phase interfaces, and (Hi) the macroscopic evolution equations in the limit e — > 0. 

To keep the presentation as simple as possible, we use a finite dimensional lattice with j = 
—N,...,N and close the resulting ODE system by imposing homogeneous Neumann boundary con- 
ditions 

u- N -i(t) = u- N (t), u N +i(t) = u N (t). (2.1) 



The natural scaling parameter on such a finite lattice is e = 1/N, that means the macroscopic 
space variable £ takes values in the interval [—1,1]. The numerical simulations presented below are 
computed by the explicit Euler scheme, where the time step size is chosen sufficiently small so that 
the energy, see (2.2) below, is strictly decreasing. Moreover, all simulations are performed with 

2(1 -u 2 ) 2 16« 



1 + n 2 ' v ; (1 + u 2 ) 2 ' 

which is convenient for numerical computations since the linear growth of <£' for \u\ — > oo allows to 
use a relatively large time-step size. 

2.1 Gradient flow structure and onset of macroscopic data 

The lattice diffusion (1.1) can be regarded as a discrete analogue to the H~ ^gradient flow of a 
nonlinear bulk energy. More precisely, defining the energy 

N 

8{u) := e Y, *(«j) ( 2 - 2 ) 

j=-N 



and the metric potential 

N 

2 



'R-i'u) = - 2_, (^+ v j) with — Avj = iij for j = —N, . . . , N and v±(jv+i) = V ±N 



we readily verify - using discrete integration by parts along with the boundary conditions (2.1) - 
that (1.1) is equivalent to 

duK(u) + d u £(u) = 0, 

where the metric tensor dyTZ is formally given by (— A)~ . In particular, we obtain the energy 
balance 

f = e*f T = s 2 V, V{u) := e^ £ (V + #'K)) 2 , 

j=-N 

where the dissipation T> gives the squared and rescaled length of the energy gradient with respect 
to the metric induced by 1Z. Notice that £ , 1Z, and T> are scaled macroscopically, that means the 
identification (1.4) implies 

£(U) = j HU)dt V(U) = j (%$'(C/)) 2 d£ (2.3) 

as well as 

K(d T U) = l (dt:V) 2 dt with -dlV = d T U and d(V\ i=±1 = (2.4) 

provided that U is sufficiently smooth with respect to r and £. The formal gradient flow corresponding 
to (2.3) and (2.4), however, is the ill-posed PDE (1.6) and can hence not govern the limit dynamics. 
A further important observation is that the nonlinear lattice (1.1) admits a comparison principle 
on the increasing branches of <£>'. Specifically, using standard arguments for ODEs we easily show 

sup sup Uj(t) < max{« , sup Uj(0)\, inf inf Uj(t) > min{u#, inf Uj(0)\. 

t>0\j\<N " \j\<N *>0|j|<Af \j\<N 
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Figure 2.1: Illustration of the initial transient regime for N = 50, using snapshots of Uj and pj (left and 
n^/ii column, respectively) against £ = ej at time r = (gray) and a short time r = 0+ afterwards (black). 



The shadowed regions in the left column indicate the intervals 



and 



u # ,u 



#1; the dotted horizontal 



lines in the right column represent the critical values {p*,p*}. The vertical lines describe the macroscopic 
phase interfaces. Top row. In this example, the initial data Uj{0) at r = do not penetrate the spinodal 
interval. The lattice data Uj and pj at r — 0+ therefore resemble macroscopic functions U and P, which are 
discontinuous and continuous, respectively, at the phase interface. Bottom row. Part of the initial data Uj(0) 
are now taken from the spinodal interval. At time r = 0+, the lattice data Uj can no longer be described 
by macroscopic functions but only by Young measures. The data pj, however, still approximate a continuous 
macroscopic function P. 



This implies £{t) = O(eN) = 0(1) for all t > and hence f™ V(u(t)) dt 
initial data Uj(Q) are bounded independently of e. 



O(l) provided that the 



We next discuss the small-time dynamics. For < e <C 1 the initial evolution of (1.1) is related to a 
very fast transient regime - where 'fast' refers to the macroscopic time r - during which the system 
quickly approaches a state with macroscopic behavior of pj and dissipation of order O(l). This is 
illustrated in Figure 2.1, which depicts two prototypical examples of initial data at r = along with 
the state of the system at a small macroscopic time r = 0+ afterwards. 

In the top row of Figure 2.1, we start with microscopic (i.e., oscillatory) initial data that are 
confined to the two stable regions u G (— oo, u*) and u G (u* , +oo). Due to the oscillations, the initial 
dissipation is of order 0(e^ 1 N) = 0(e -2 ). The initial energy gradient is therefore also very large 
and drives the system rapidly. At time r = 0+ ~ 0(e 2 ), the dissipation and the energy gradient 
eventually become of order O(l). This implies that the discrete data pj resemble a macroscopic 
function P that admits a weak spatial derivative and is hence continuous with respect to £. We also 
observe that the discrete data Uj at r = 0+ approximate a piecewise continuous function U which 
satisfies P = $'([/) and jumps across the interface located at £ = 0. In particular, the phase fraction 
H defined by 

/« := X[u*,+oo){U) - X(-oo,u t )(U) 

takes the values —1 and +1 outside of the interface, where \J denotes the indicator function of the 
interval J. 

The second example, see the bottom row in Figure 2.1, is different since now some of the initial 
data Uj(0) belong to the spinodal interval [u*, u*], in which the discrete diffusion coefficient is negative. 
In the numerical simulation we therefore observe that each of those Uj quickly leaves the spinodal 
interval (spinodal decomposition) and that the data for adjacent j can be attracted by different 
stable regions. In particular, the data at time r = 0+ ~ 0(e 2 ) exhibit a phase interface near 
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£ ~ —0.23, in the sense that Uj is non-oscillatory on the left but highly oscillatory on the right of 
the interface. The fine structure of these oscillations depends on the microscopic details and each 
reasonable macroscopic theory must describe them in terms of a Young measure v = v(t, £, du), 
which provides a probability distribution with respect to u for any macroscopic point (r, £). The 
discrete data pj, however, still resemble a macroscopic function P because otherwise the dissipation 
could not be of order O(l). Since pj and Uj are coupled by <&', we then conclude that the u-support 
of the Young measure v consists of only two points. This reads 

i/(r, £, du) = <5/3_(p(t,0)( dw ) + g 5 /MP(U))( du )- ( 2 - 5 ) 

Here, <5g(du) is the Dirac distribution at /3, the functions /3_, /3+ denote the two stable branches of 
the inverse of <J>', and the phase fraction \x takes values in [— 1, +1]. 

The simulations from Figure 2.1 reveal that there exist (at least) two different types of macroscopic 
phase interfaces: Type-I interfaces separate regions where the microscopic data Uj are confined to 
either one of the phases (— oo,u*) and (u*,+oo), whereas type-II interfaces describe that the lattice 
data oscillate between the two phases on at least one side of the interface. Below we argue that type-I 
interfaces can be described by a free boundary value problems, which exhibits hysteresis and involves 
only the macroscopic fields P and \i G {—1, +!}• Type-II interfaces, however, are more complicated 
and their investigation is postponed to future research. 

2.2 Examples of macroscopic interface dynamics 

In this section we study the dynamics of type-I interfaces in numerical simulations. In particular, we 
investigate the macroscopic jump conditions across such interfaces and provide numerical evidence 
for pinning, depinning, and annihilation. At the end we also present an example of a type-II interface. 
For simplicity, and in view of the discussion in the previous section, we always impose initial data 
Uj(0) such that Pj(0) = $>'(uj(0)) resemble a macroscopic function £ t- >■ P(0, £). In all simulations 
we observe - for, loosely speaking, most of the macroscopic times r > - that the discrete data 
Pj (t/g 2 ) approximate a macroscopic function £ t- > P(t, £). We therefore expect that the macroscopic 
limit e — > can in fact be characterized by a PDE for P and the phase field fi, or equivalently, in 
terms of a free boundary problem for P = $>'(U) and the interface curves. There exist, however, 
small macroscopic times intervals in which the discrete data pj(t) exhibit strong temporal and spatial 
fluctuations near a moving phase interface. These fluctuations are discussed in the next section. 

Figure 2.2 depicts lattice simulations with two different values of TV, showing that the macroscopic 
plots of the discrete data are basically independent of TV. In this example, we initialize two macro- 
scopic phase interfaces which are located at £ = and £ = 0.6 and separate regions with U > vft, 
U G [?/#,«*], and U G [u*,u^]. The first interface moves to the right while the second one clearly 
keeps its initial position. At the later time r ~ 0.18 both interfaces annihilate each other in a colli- 
sion process, see also Figure 2.4, and the macroscopic evolution afterwards is governed by nonlinear 
diffusion inside the phase U G [u*,oo). Figure 2.3 provides numerical evidence for the jump rules 
across the interface: The moving interface is driven by a jump in d^P whereas P is smooth across the 
standing interface. Moreover, while P = p* holds on the moving interface, P evolves on the standing 
interface and takes values in [p*,p*]. 

A further dynamical feature of the lattice (1.1), namely the pinning of interfaces, is illustrated 
in Figure 2.5. At time r = 0, we initialize a single macroscopic interface that separates regions with 
U > u* and U G [u#, u*], where the data are chosen such that P > p* and P < p* holds locally on 
the left and on the right of the interface, respectively. This interface starts propagating to the right 
but stops moving at r ~ 0.02 because the bulk diffusion behind the interfaces enforces P < p* for 
r > 0.02. The inverse process, that is the depinning of interfaces, is shown in Figure 2.6. There, 
a single macroscopic interface is initially at rest with P G (p*,p*) but propagates with P = p* for 
r > 0.03. 
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Figure 2.2: Annihilation of a moving and a standing type-I interface for two different simulations with 
TV = 200 (thick gray curves) and N = 500 (thin black curves). The snapshots show Uj against the macroscopic 
position £ = ej at fixed macroscopic times r; the vertical lines indicate the interface positions and the two 
shaded regions represent the intervals [u*,u*] and [u#,u*]. 
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Figure 2.3: Snapshots of Pj for the example from Figure 2.2. The horizontal line represents p = p* 
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Figure 2.4: Macroscopic positions of the interfaces and the energy as function of r for the example from 
Figure 2.2. 



The last example is shown in Figure 2.7 and concerns the dynamics of type-II interfaces. The 
discrete initial data Uj(0) resemble a smooth function U that penetrates the spinodal region, and the 
initial transient regime therefore leads to spinodal decomposition. This means the lattice dynamics 
forms a phase interface that separates a region with U > u* from a region with strong microscopic 
oscillations, where the latter can be regarded as an approximation of a nontrivial Young measure 
v. The interface is initially at rest but depins at r ~ 0.02 and propagates into the oscillatory phase 
afterwards. The dynamics of type-II interfaces are not yet well-understood. In particular, since 
there exist many microscopic realizations of a given Young measure v, it is not clear whether the 
macroscopic evolution can be completely characterized in terms of the fields P and (i G [—!,+!], 
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Figure 2.5: Pinning of a type-I interface with N = 400. The vertical lines indicate the initial and the current 
position of the interface. 
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Figure 2.6: Dcpinning of a type-I interface with N — 500. 

or whether further macroscopic quantities are needed. All subsequent considerations are therefore 
restricted to type-I interfaces. 

2.3 Microscopic dynamics of phase interfaces 

To conclude our numerical investigations, we now discuss the microscopic dynamics of moving type- 
I interfaces in greater detail. In this way we not only obtain a better understanding of the lattice 
dynamics but also gain some inside into the analytical problems that must be addressed when passing 
to the limit e — > 0. 

For a moving type-I interface that is isolated - i. e., sufficiently far from any other interface and 
the boundary of the computational domain - the key numerical observations are illustrated in Figures 
2.8 and 2.9 and can be summarized as follows. 

1. Sequentiality. At each time there exists at most one index k such that Uk is inside the spinodal 
interval [u*,u*]. In other words, the interface moves because the Uj's undergo the phase tran- 
sition one after another, where 'phase transition' means passage through the spinodal interval. 
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Figure 2.7: Type-II interface with N = 50 that propagates into a region with microscopic oscillations. In 
the macroscopic limit, these oscillations cannot be described by a function U but only in terms of a Young 
measure u, see 2.5. 

2. Fluctuations near the interface. Each phase transition produces strong microscopic fluctuations. 
These fluctuations are initially very localized, correspond to relatively large dissipation, and 
are spread out before the next phase transition occurs. 

3. Temporal scale separation. The time needed to spread out the fluctuations is as least as large 
as the time to undergo a phase transition. Moreover, both times are typically much smaller 
than the time between two subsequent phase transition. 

The fundamental sequentiality property can - at least in an idealized single-interface setting - be 
derived from elementary comparison principles for ODEs. The underlying idea is that as long as Uk 
is inside of [it*, it*], the lattice equation ensures that any other Uj cannot enter the spinodal interval. 
For the piecewise quadratic potential, we employ a similar argument in the proof of Theorem 3.2 in 
order to show the persistence of single-interface data. 

The microscopic fluctuations are much harder to describe rigorously. More precisely, although it 
is relatively simple to understand the onset of fluctuations heuristically, it is not obvious, at least to 
the authors, how to estimate their spatial and temporal decay in the case of a generic double-well 
potential <&. For the piecewise quadratic potential, however, the fluctuations can be controlled by 
splitting p = $'(u) into a regular part related to linear diffusion and a sum over delayed and shifted 
variants of the discrete heat kernel, see the discussion in §3. 

A further challenge for any rigorous treatment is to give a suitable description of the different time 
scales in the problem. For instance, in order to guarantee that each interface propagates with finite 
speed on the macroscopic scale, one has to show that the microscopic time between two adjacent 
phase transitions is of order 0(e~ 1 ). Moreover, proving that any limit function P is in fact continuous 
at the interface requires to show that both the time for each phase transition and the decay time 
of the fluctuations are much smaller than 0{e~ 1 ). For piecewise quadratic $, the problem is again 
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Figure 2.8: Microscopic dynamics of a moving typc-I interface: snapshots of Uj(t) (black) and pj (gray, 
rcscalcd) against j G Z at six non-equidistant times; the vertical line indicates the interface position at j = k. 
At time t\, the value Uk has just crossed u* from below. Afterwards, it passes rapidly through the spinodal 
interval [u*, u*] and evokes microscopic fluctuations which are of order 1 but localized near j rj k. At time £3, 
the value Uk leaves the spinodal interval but the fluctuations are still present and not spread out before the 
time £5. Between £5 and te, the data change only little but prepare u^+i for the next phase transition, that 
means Uk{te) > u* and Uk+i{te) = u*. 




Figure 2.9: Microscopic dynamics of a type-I interface that moves from j = 1 to j = 8. The black curves 
represent the time trajectory of U4, and the gray curves describe the (rescaled) dissipation of the system. Top. 
$ is a smooth double well-potential and each interval i| contains all times at which Uk is inside of the spinodal 
interval. Bottom. For the piecewise quadratic potential, each interval 1%. is degenerate and consists of a single 
time t* k with Uk{t* k ) = 0. 

much simpler. At first, phase transitions take place at precise times due to the degenerate spinodal 
region, and second, all other time scales can be related to the properties of the discrete heat kernel. 

2.4 Effective evolution equations for the macroscopic dynamics 

We now derive the free boundary problem for the dynamics of type-I interfaces on a heuristic level. 
Since our arguments are very similar to those for the viscous approximation, we only sketch the main 
ideas and refer to [EP04] for more details. 

We first suppose that the lattice data Uj and pj converge strongly as e — > to macroscopic 
functions U and P which are sufficiently regular. For simplicity we also assume that there is only a 
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single interface located at £*(t) and that U satisfies the phase condition 

U>u* for £ < £* an U < u* for £ > £*. (2.6) 

Under these assumptions, and using the weak formulation of (1.5), we readily verify that that macro- 
scopic evolution is governed by bulk diffusion via 

d T U = dfP with P = &(U) for all r > 0, £ ^ f*(r) (2.7) 

and the generalized Stefan condition 

\[P]\=0 and ^|[t/] + |[%P]|=0 for £ = £*(r), (2.8) 

which ensures that (2.7) holds in a distributional sense even across of the interface. Here \[U} denotes 
as usual the jump of U across the interface, this means 

[17] (r) = U + (t) - C7_(r), U ± (r) := lim [7(r, f (r) ± h). 

The nontrivial part is to identify a further dynamical interface condition since (2.7) and (2.8) do 
not the determine the evolution of £* completely. In view of the numerical results, we propose the 
following hysteretic flow rule, see Figures 1.2 and 1.3: At almost each time r > 0, the interface is 
either 

1. standing with \[P}\ = |[<%P| = -^C = and P e [p*,p*], or 

2. propagating into U < u* with P = p*, -g^£* > 0, and \[U]\ = u* — u* < 0, or 

3. propagating into U > u* with P = p*, -g^£* < 0, and [C7| = u# — u* < 0. 

Notice that the combination of phase condition, bulk diffusion, Stefan condition, and flow rule pro- 
vides - at least on a formal level - a well-posed free boundary value problem and that the limit model 
can easily be generalized to an interface with U < u* for x < £* and U > u* for £ > £*, and to the 
case of finitely many phase interfaces. 

The interface laws can also be derived in a more sophisticated way. For the viscous approximation, 
it has been shown in [EP04] - also assuming sufficiently strong convergence as e — > - that any limit 
solution must satisfy the distributional entropy condition 

d T ^{U) - %(T(P)%P) < (2.9) 

where T is an arbitrary but increasing function and \& defined by \£'([7) = T(&(U)). This family 
of local entropy inequalities can also be established in the lattice case. In fact, the scaled equation 
(1.5) implies 

d T *{U) = T(P) A £ P = V +£ (T(P)V_ e p) - (V +£ T(P)) (V +£ P) < V +£ (t(P)V_ £ P 

due to T' > 0, and passing to e — > in the weak formulation of (1.5) we get (2.9). The key 
observation is that (2.9) combined with (2.8) yields 



[*(C7)]| - T(P)\[U]\ < and hence ^- ( f + T($'(s)) - T(P)ds J < 0, 



where c7_ > [7+ and P_ = P+ = P hold thanks to (2.6) and (2.8). The flow rule now follows from 
evaluating the latter inequality for appropriately chosen functions T, see again [EP04] for the details. 

In the special case that $ is given by the piecewise quadratic potential (1.11), the bulk diffusion can 
be written as 

d T P = dfP for all r>0, £^£*(r), 
13 



whereas the Stefan condition simplifies to 

\[P}\ = and 2^ = \[d^P}\ for £ = £*(r). 

In the space of distributions, both equations can be combined to 

a r (P + ft) = d}P for all r > 0, £ G R, (2.10) 

where the phase field 

M(r,6:=sgn(r(r)-e) (2.11) 

is well-defined outside of the interface and takes values in { — 1, +1}- We mention that the hysteretic 
flow rule can be encoded as 

fi = M[P] (2.12) 

where M is a particular relay operator, see for instance [Vis94, BS96], and that well-posedness of the 
initial value problem for (2.10) and (2.12) has been proven in [Vis06]. For our purposes, however, 
the differential relations from above are more convenient than the discontinuous hysteresis operator 
A4. In particular, under the sharpened phase condition 

U(t, > u # = -2 for £>T(t), 

which implies that the interface is either at rest or propagates to the right, the flow rule is equivalent 
to 

P(t,C (r)) € [-1,+1] with P(t, e*(r))=+l for -^(r) > 0. (2.13) 

In §3 we prove that the simplified free boundary value problem (2.10) with (2.11) and (2.13) indeed 
governs the parabolic scaling limit of (1.1) with (1.11) provided that we impose macroscopic single- 
interface initial data as described in Assumption 3.6. 

3 Rigorous analysis for the piecewise quadratic potential 

In this section we characterize the lattice dynamics with piecewise quadratic potential. In particular, 
from now on we suppose that 

$>'(u) = u — sgnu, (3-1) 

where we define the sign function by 

/ -1 for u < 0, . . 

SgnU= ( +1 for n>0. (3 ' 2) 

The condition sgnO = +1 is essential for establishing global existence and uniqueness of single- 
interface solutions, see the proof of Theorem 3.2 and the remark afterwards. 

Our goal in this section is to prove (i) that a certain class of well-prepared microscopic initial data 
give rise to a single phase interface that moves in a certain direction, and (ii) that the macroscopic 
dynamics for e — > can in fact be described by the simplified free boundary problem from §2.4. 
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3.1 Existence of single-interface solutions 

In order to make precise what we mean by single-interface solution, we define for each k G Z the 
state space 



Xk 



u G 



< inf Uj < sup Uj < oo, — 2 < inf Uj < sup Uj < > 
3< k j<k i> k j>k J 



see Figure 3.1 for an illustration. As shown below, the dynamics for initial data chosen from Xk is as 
follows: The system stays inside Xk until some maximal time tt > with linnet* Uk(t) = at which 
Uk undergoes a phase transition by crossing the spinodal value from below. For times t > t k , the 
system evolves inside of Xk+i until Uk+\ changes from negative to positive phase at time t* k+l > t* k . 
By iteration we therefore find a phase interface that moves to the right, where we allow for tt 
to account for standing interfaces. 



DC 
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Figure 3.1: States from Xk have a single phase interface at k (left panel, Uk-i and Uk arc highlighted). Under 
the dynamics, this interface is either stationary for all times or shifts eventually to the right due to a phase 
transition of Uk, because the system can reach the boundary dXk only via Uk = 0. When this happens at time 
t* k (right panel), we also have uu~i(t* k ) > 2 and ilk jumps from Uk(t* k — 0) > to Uk(t* k + 0) > 4. Afterwards 
the state of the system belongs to Xk+i until Uk+i undergoes a phase transition at time t* k+1 . 

Definition 3.1 (single-interface solution). A continuous function u: [0, oo) —> ^°°(Z) is called single- 
interface solution to (1.1) with potential (3.1) if there exists k\ £ Z along with a sequence (t* k ) k>k C 
(0, oo] such that the following conditions are satisfied with t* k _ 1 := 0: 

1. For all k > k\, we have either t* k _ l = oo or t k _ 1 < t* k , 

2. The function u solves iij = A &(uj) for all times t G [0, oo) \ {t k : k > ki}, 

3. We have Uk(t* k ) = and linn,^* Uk(t) > for all k > k\ with t* k < oo, 

4. The function u takes values in Xk on each time interval {t*k_i,t* k ) with t k _-y < oo. 

Moreover, a single-interface solution with t* k = oo for all k > k\ is called a standing-interface solution. 

For single- interface solutions in the sense of Definition 3.1, any Uk is continuous at the phase 
transition time t k , whereas pk jumps at t* k from +1 to —1. Notice that it is the other way around 
in the macroscopic limit: when e tends to 0, p becomes continuous with respect to both r and £, 
whereas u jumps from to +2 at the phase interface, see Figure 2.9 for an illustration. 

As first main result we prove that single- interface solutions exist. 

Theorem 3.2 (existence and uniqueness of single-interface solutions). For given initial data in Xk i; 
there exists a unique single-interface solution to (1.1) with (3.1). Moreover, this solution satisfies 

sup supuj(t) < D := max {2, supiij(O)}, 

te[o,oo) jez jez 
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and we have 



u k -M) > 2, 4 +1 - t* k > In Jyr^l, «*(**) = Jim «*(*) = 4 + lim ti fc (t) > 4 

for all k>k\ with tt < oo. 

Proof. We assume, without loss of generality, that k\ = 1 and derive all claims by induction with 
respect to k > fci. 

Dynamics inside of X\: For states u G Xi we have sgnu.,- = — sgn (j — 1) and hence 

iij = A &(uj) = A Uj + 25° - 2(5), (3.3) 

where <5^ is the Kronecker delta. In particular, the right hand side of (3.3) is Lipschitz continuous 
with respect to the £°°-norm of u and this implies the local existence of a unique solution that is 
smooth in time and takes values in the open set X\. From the definition of X\ we also infer that 

< Uj(t) + 2uj{t) < 2D iorj<l, 

< Uj(t) + 2uj(t) < D + 2 for .7=0, 

-4 < Uj(t) + 2uj(t) for j = l, 

-4 < Uj(t) + 2 Uj {t) < forj>l, 

and the comparison principle for scalar ODEs gives 

e~ 2t Uj(0) < Uj(t) < e^ 2t Uj(0) + (l-e~ 2t )D for j < 1, 

e- 2t u j (0)-2(l-e- 2t ) < Uj (t) for j > 1, 

Uj(t) < e- 2t Uj(0) (oij>l. 



Using these estimates, we now conclude that the local solution with values in X\ either exists for all 
times (in which case we set t\ := oo) or reaches the boundary of X\ at some time < t\ < oo via 
ui(t\) = 0. Moreover, since both u\ and iii are continuous at any t G [0, t\) we have \m\ t/ A t* u\{t) > 0. 
Phase transition at t\: Now suppose that t\ < oo. From (3.3) and the above estimate for Uj(t) 
we conclude that Uj{t) converges for any j G Z as t /* t\ to some well-defined limit Uj(tl), where 
ui(t\) = as well as 

< inf Uj{tl) < supuj(tl) < D, -2 < mfuj(tl) < supuj(tl) < 0. 

We therefore have 

sup lAS'^- (£*))| <C 

for some constant C. On the other hand, uo(i*) > and —2 < ^2(^1) < ensure that 
A$'(m(^)) = (uo(tl) + u 2 (tl)) - (sgnu ( y t* 1 ) + sgnu2(t* 1 )-2sgn0) >0 



thanks to sgnO = +1. These results reveal that the vector field of the dynamical system u = A <&'(u) 
is well defined in the state «(£*) G dX\ (18X2 and points into Xi, where dX^ denotes the topological 
boundary of X k in £°°(Z). We therefore find a time t* > t\ such that the local solution to the 
microscopic dynamics from the first step has unique continuation to the time interval [0,t*) with 
u(t) G X2 for all t G (£*,£*). By construction, Uj and iij are continuous at t\ for j / 1, whereas 

ui(t\) = lim ui(t) = 0, uiitV) = lim iiUt) = 4 + lim uUt) > 4. 

thanks to sgnui(t^) = lim^i* sgnui(t) = +1 and lim^* ^i(*i) = —1- Finally, repeating all ar- 
guments from the first step we show that this solution in X2 exists uniquely until a time t\ with 
t\ < t\ < 00 and U2(^2) = in case of t^ < 00. 
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Condition for uo(i*) and lower bound for t* 2 — t\: For < t < t* we infer from (3.3) that 

«i(t) = (u (t) + u 2 (t) - 2«i(t)) - ((+1) + (-1) - 2(-l)) < uo(t) + u 2 (t) - 2«i(t) - 2. 

Passing to the limit £ /* t* we therefore get 

< ui(tj) < «o(<i) + «2(*i) - 2 < wo(*i) - 2 

and hence uo(t*) > 2. In the same way we prove u\{t^) > 2 and conclude that the interface can shift 
from k = 1 to k = 2 at time t\ only if 

2< wi(*5). 

Moreover, u(t) G X 2 for all t G (i*,^) implies 

«i(t) = («o(*) + «2(*) " 2wi(t)) - ((+1) + (-1) - 2(+l)) < D + 2 - 2ui(t). 
By the comparison principle we therefore have 

t .\D + 2 



ni(t)<(l-e- 2<+2 ^) 



and setting t = t\ we obtain the estimate for ££ — **• 

Final step: All arguments derived above can easily be iterated, so the assertions of Theorem 3.2 
are proved by induction. In particular, the lower bound for t* k+1 — t* k ensures that the solution exists 
for all times t > 0. □ 

Remark 3.3. Theorem 3.2 can be adapted to finite lattice systems with j = 1, . . . , M provided that 
the corresponding ODE system is closed by imposing 

1. either homogeneous Neumann boundary conditions via uo(t) = u±(t) and UM+i{t) = UM(t), 

2. or inhomogeneous Dirichlet boundary conditions via uo(t) = c\ > and UM+i(t) = c 2 < 0, 

where j = and j = M + 1 are the ghost indices. 

The proof of Theorem 3.2 reveals that the condition sgnO = 1 (or equivalently, the right-sided 
continuity of <&') is crucial for the microscopic dynamics to be well defined after the phase transition 
time t* k as it guarantees that Uk{t* k ) = implies the strict inequality Uk(t k + 0) > 0. Our convention 
(3.2) is therefore implicitly but intimately related to phase interfaces that propagate into the phase 
u < (to observe phase propagation into to the other phase u > 0, one has to set sgnO = — 1). This 
subtle direction selection reflects that the spinodal region is actually an isolated point and can be 
regarded as the degenerate analogue to the hysteresis flow rule (1.10). 

Corollary 3.4 (criterion for standing interfaces). For single-interface initial data 

u(0) £ Xfcj with sup'Uj(O) < 2 

we find t* k = oo and hence a standing-interface solution. 

We finally characterize the evolution of Pj(t) = $>'(uj(t)) for single-interface solutions. The key 
observations are as follows: 

1. p solves the linear discrete heat equation p = Ap on each time interval (^'^fe+i) because 
u(t) £ Xj~ implies Pj(t) = Uj(t) + sgn (j — k) and hence Pj(t) = Uj{t). 

2. At each phase transition time tt < oo, the value of pk jumps down according to 

lim p k (t) = +1, lim p k (t) = -1 

t/t\ t\t* k 

but we have lim^t* Pj(t) = lini£\t* Pj(t) for all j ^ k. 
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The dynamics of p can therefore be interpreted as solving the linear discrete heat equation for t > 
where at each phase transition time ££ we perturb the system by adding —26^ to the current state. 
Combining this interpretation with the superposition principle we conclude that p consists of a regular 
part q and a singular part r, which account for effects of the initial data and the perturbations induced 
by the phase transitions, respectively. In particular, denoting by g the discrete heat kernel, that is 



gj (t)=--A 9j (t), 

9j(0) 



<5 U 



(3.4) 



we arrive at the following result, see also Figure 3.2. 

Corollary 3.5 (decomposition of p). For each single-interface solution we have 



with 



Pji*) = Qj( t ) + r j(. t )^ 



r i(*) := ~ 2 J2 X{t>tj}(*)i&-fc(* ~ **)' 



k>k! 



where q is the solution of the discrete heat equation with initial data qj(0) = Pj(0) and \J denotes 
the indicator function of the set J. 



n9j{t) 



29j-k(t-tl) 




Figure 3.2: Left. Cartoon of the discrete heat kernel g (the thin lines represent interpolating splines and are 
drawn for better illustration). Right. The lattice function r is at any time t given by a finite sum of shifted 
and delayed version of the discrete heat kernel. 

The decomposition formula from Corollary 3.5 is crucial for our analysis. In particular, the decay 
and continuity properties of the discrete heat kernel, see appendix A, enable us to estimate the 
impact of the sequence of singularities induced by the phase transitions. Notice however, that the 
evolution of p is still nonlinear as the phase transition times t* k depend on p itself. 



3.2 Upper bounds for the macroscopic interface speed 

In Theorem 3.2 we have shown that single-interface solutions exist. In order to pass to the macro- 
scopic limit, we must however guarantee that the corresponding macroscopic interface speed is at 
most of order Oil). More precisely, for a given macroscopic distance 8 > we have to ensure that 
the macroscopic time £ 2 t%i < is bounded from below by some constant independent of e, where [^J 
is shorthand for the integer part of x. The derivation of this lower bound is actually at the heart of 
our convergence proof and requires a careful quantitative analysis of the lattice dynamics. In this 
paper we restrict ourselves to the most simple case and suppose that the initial data are compatible 
with the macroscopic limit model. We also assume without loss of generality that the macroscopic 
interface is initially located at £ = 0. 

Assumption 3.6 (macroscopic single-interface initial data). The initial data u(0) belong to X\ and 
p(0) = $'(u(0)) satisfies 



|Apo(0)| < /?£, 
for two constants a and j3 independent of e 



sup|V + pj(0)| < as, sup | A pj (0)| < ae 
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Assumption 3.6 is motivated by the limit dynamics, see also Figure 3.3 for an illustration. In 
fact, the prototypical example for macroscopic single-interface data is 

Pj(t) = flni(ej) + ce, uj(0) = Pj(t) + sgn (-j) (3.5) 

where c £ is a constant and P; n i some given macroscopic function independent of e such that: 

1. -Pi n i is bounded and continuous for all £ G R with 

P ini (£) > 1 for £ < and - 1< P ini (£) < 1 for £ > 
and 

liminf P ini (£) > 1 and - 1 < liminf P in i(£) < limsupP in i(0 < + 1 - 

£-►-00 C^ + oo ' ^+oo 

2. Pi n i is twice continuously differentiate for £ < and £ > with 

:= lim |i^(+0 - PU-0\ <™, ol := sup |^ ni (£)| + |i&(£)| < °o, 

3. c e > is sufficiently small with c e — >■ as e — > and chosen such that t\ > 0, that means 
Uj(0) > for all j < and Uj(0) < for all j > 1. 

Corollary 3.7 (bounds for the regular part g). VFe /iave 

|V + ^(t)| < ae and \qj(t)\ = \Aqj(t)\ < ae 2 + (3eg {t) 

for all t > 0, j G Z, and e > 0. 

Proof. By construction and Assumption 3.6 we have 

£<&■(*) = A^-(t), |^(0)| = 1^(0)1 = |^(0)| = |A Pj (0)| < ae 2 + (3e5° 

as well as 

f t V +qj (t) = AV+ qj (t), |V+^(0)| = |V+Pj(0)| < ae. 

The claim now follows using both the superposition and the maximum principle for the discrete heat 
equation. □ 

We remark that the assertions of Corollary 3.7 are sufficient for showing that the macroscopic 
interface speed is bounded. All results derived below therefore hold (with different constants) even 
in the case that 

1. \V+Pj{0)\ < ae and |Apj(0)| < pe for all j G Z, 

2. \Apj(0)\< ae 2 for almost all j G Z, 

that means the derivative of the function Pi n i from (3.5) can even be discontinuous at finitely many 
points. 

We are now able to derive upper bounds for the macroscopic interface speed. To this end we prove that 
Assumption 3.6 implies that the time t* k+1 — t* k between two adjacent phase transitions is mesoscopic 
as it can be bounded from below by 1/e (recall that Theorem 3.2 provides only t* k+l — t* k > 2/(D — 2)). 
Our considerations are based on the estimate 

2<p k (t* k+1 )-p k (t* k )= l + p k {t)dt= I ' q k (t)dt-2j2 l + 9k-n{t-t* n )dt, (3.6) 

^t ^t n=l ^ l l 

which follows from combining the conditions Uk(t* k ) = and u k (t k+1 ) > 2 with the representation 
formula from Corollary 3.5. 
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Figure 3.3: For initial data as in Assumption 3.6, the discrete data Pj(0) resemble a macroscopic function 
Pi n i that is continuous and piecewise twice continuously diffcrentiable. In particular, the initial velocities 
Uj(0) = Pj{0) satisfy po(0) = fie and |pj(0)| < ae 2 for j ^ 0. 

Lemma 3.8 (refined lower bound for the time between two phase transitions). For any Tfi n > 
there is a constant d* > 0, which depends on a, 6, and Tfi n , along with a constant < e* < 1, which 
depends only on a and 8, such that 

e{t* k+1 -tl) >2eU 

holds for all k > 1 with < tl < Tg n /e 2 provided that e < e*. 

Proof. In what follows we consider k > 1 with t* k < Tfi n /e 2 and denote by c and C generic constants 
independent of a, 8, and e. 

Weaker variant of (3.6): In order to study the implications of (3.6), we first simplify the right 
hand side as follows. By Lemma A. 2 we have 



fc+i 



Y, I 9k-n{t-t* n )dt< 

whereas Corollary 3.7 provides 



Yl I 9o(t-t* n )dt = ^2go(t* k -t* n )-^2g (t* k+1 -t* n ). 

n=l *fc n=l n=l 



'fe+l 



g fc (t)dt< 



fc+i 



ae 



+ 8eg (t)) dt 



Rearranging terms and using go(0) = 1, inequality (3.6) becomes 

k ., ,.+* fe-1 



!>(**+!-£) ^ £ //^ {^ 2 + Beg Q (t))dt + Y,9o(tt-tl), 

n=l •**% n=l 



and writing s k '■= t k+1 — t* k we arrive at 



2^5o(sfe H \-s n ) < - I (ae 2 + 3eg (t)) dt + 2_,#o(sfc-i H U„ 

n=l ^*fc n=l 

In the remainder of this proof we transform this inequality into a lower bound for s k - 
Estimate for k = 1: Inequality (3.7) combined with Lemma A.l provides 



(3.7) 



9o(si) 



< \ / (ae 2 + 8eg (t)) dt < \ / (ae 2 + (3eg (t)) dt 

Jt\ Jo 



and we deduce that there exists e* > 0, which depends on a and 8, such that s\ > 1 for all e < e*. 
We employ Lemma A. 2 again to estimate 



and this gives 



<C(ae 2 + /3eVsT) 



esi > d\ > 
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for some constant d\ which depends on a and /3 but not on e. For /3 = we even find e 4 ' 3 si > d\. 
Estimate for k > 1 in a special case: For the following considerations we further suppose that 

s k < min s n . 

n=l,...,k— 1 

Since go is strictly decreasing, we therefore get 

go(sk H h s n+ i) > 5o(sfe-i H 1- s n ) 

for all n = 1, . . . , k — 1, and hence 

3o0fc H h si) < - / (ae 2 + /3eg (t)) dt, 

2Jt l 

thanks to (3.7). This implies fk{sk) < for 

1 f l l +s 
fk(s) := go(s + t%)-- (ae 2 + f3eg (t)) dt, 

where we used that <?o( s fc + - ■ ■ + si) = go(sk + t* k — t*) > 9o( s k + *&)■ The properties of go, see Lemma 
A.l once more, guarantee that the function /& is convex, continuous, and strictly decreasing in s. 
Therefore, and in view of /&(()) > and hin^^oo fk(s) = — oo, we conclude that fk has a unique zero 
Sfc with Sk < Sk- Due to convexity of fk we also find 

Sk < Sk < s k , 

where §k '■= —fk(0)/f' k (0) is the first approximation to s~k when starting the Newton algorithm at 
s = 0. By direct computations we now verify 

s 90(tt) 

Sk ~ 



-9o(tt) + ¥ + f9o(tt 
and using tt > si > 1 as well as Lemma A. 2 we obtain 



and hence 

es k > — -rr 2 =: d 2 

d^ + ocr^ + P 

where we used that t* k > s\ > d\/e and t* k < Tfi n /e 2 . 

Estimate for k > 1 in the general case: We have established the estimate es\ > d\ as well as the 
implication 

es k < minjesi, • • • , esfc-i} =>■ es k > d 2 , 
and the desired estimate for Sk follows with d* := ^ minjcii,^} by induction. D 

From Lemma 3.8 we immediately obtain t* k > 2kd*/e and we deduce for each macroscopic time 
Tfi n that at most T^ n /(2ed^) phase transitions can happen for r < Tfi n , shifting the interface over a 
macroscopic distance smaller than Tfi n /(2d,,). 

We conclude this section with some comments concerning the microscopic fluctuations caused by 
the phase transition of Uk at tt. The properties of the discrete heat kernel imply that the amplitude 
as well the inverse length of the effective spatial support scale with 



£ 2 + ( T _ e%) 



21 



which decays quite rapidly within a macroscopic time of order e 2 but much slower afterwards. In 
particular, when Uk+i undergoes the next phase transition at time t* k , 1; the fluctuations evoked by 
Uk have reached an amplitude of order e 1 ' 2 and spread over a macroscopic length of order e^ 1 ' 2 . 
Similarly, the amplitude of the velocity fluctuations at time ££ +1 is of order e 3 ' 2 . 

These scaling arguments, especially the fractional powers of e, reveal that macroscopic single- 
interface data are not invariant under the dynamics. In other words, the lattice data at times t^,t* k 
- this means just before the phase transitions - do not satisfy Assumption 3.6, and we conclude that 
interface propagation in discrete forward-backward diffusion equations is not only a two-scale but a 
genuine multi-scale problem. 

3.3 Macroscopic continuity and compactness results 

In order to pass to the macroscopic limit e — > 0, we regard the discrete data Pj(t), qj(t), and Tj(t) as 
piecewise constant functions with respect to the macroscopic coordinates (r, £). More precisely, we 

set 

Pe (r, i + C) := VHe (r/e 2 ) for all r > 0, £ G eZ, C G [-e/2, +e/2) 
and define Qe and i? £ by analogous formulas. We further introduce the macroscopic interface position 



Ce(r):=eJ2kX[tt,t t+1 )(r/e 2 ) for all r > 0, 



fe=i 

as a piecewise constant function in time that jumps at the macroscopic phase transitions times defined 
by 

T k,e := £2 t*k f° r au k G Z. 

In what follows we fix a macroscopic final time < rg n < oo and wish to pass to the limit e — > on 
the macroscopic time-space domain 

Q:=IxR, 7:=[0,rfin]. 

To this end, we recall that Lemma 3.8 provides constants d* > and < e* < 1 such that 

A '' S 5i- Jl^.-t,>* (3-8) 

holds for all < e < e*, where 

i^ e := max fceZ : t£ e < Tfi n > 

is the number phase transitions taking place in the microscopic time interval [0, T^ n /e 2 \. Notice that 
d*, £*, and all constants derived below depend on Tfi n and on the initial data via a and j3. 

Our first result in this section concerns the compactness of the discrete interface curves £*. 

Lemma 3.9 (compactness of interface curves). The family (£*)o<e< £ * s compact in L°°(I) and eac/i 
/imii curve is Lipschitz continuous. 

Proof. We define a piecewise linear function £* on [0, rg n ] by the condition 

for all A; = 1, . . . , K £ and ^*(rfi n ) = £*( T x e ). We readily check that 

|C(r)-C(r)|<£, 0<^eW<2^ 

for almost all r G [0,Tfi n ], and conclude that the family (£g) 0< . < is bounded in C 0,1 (I). All claims 
hence follow from standard results in real analysis. □ 
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Our main goal this section is to derive compactness result for R £ and Q £ that imply (i) the 
existence of pointwise convergent subsequences, and (ii) the continuity of any limit function. In a 
preparatory step we next derive an auxiliary result for piecewise constant functions F £ on Q revealing 
that for L°°-compactness it is sufficient to establish uniform Holder estimates with respect torsi 
and discrete £ G eZ. Here piecewise constant means, as above, continuous with respect to r but 
spatially constant in each interval £ G [ej — e/2, ej + e/2), j G Z. Our auxiliary result and its proof 
are straight forward and can very likely be found somewhere in the literature on numerical analysis 
(though we are not aware of any reference) . 

Lemma 3.10 (compactness criterion for piecewise constant functions). Let (F £ ) 0<£<£ be a family of 
bounded and piecewise constant functions on f2, and suppose that there exist constants 71,72 G (0, 1] 
and C > such that 

\f £ (t 2 , 6) - F e ( Tl , 6)| < C(\t 2 - n| 71 + |& - 6| 72 ) 

holds for all < e < e* ; every t±,T2 G /, and all £1, £2 G eZ. Then this family is compact in L°°(f2) 
and any limit function is locally Holder continuous with exponent min{7i,72}. 

Proof For each e we define a piecewise linear function F £ by 

F £ (r, I) = F £ (r, £) for all r G /, £ G eZ, 

and our assumptions yield the Holder continuity of F £ with respect to time, that means 

\F £ (t 2 , 0-F £ (n,0\ <C|r 2 -ri| 71 

for all £ G R and ri,T2 G i". Moreover, since F £ is piecewise linear with respect to £ and due to 

F £ (t,Z±e)-F £ (t, I) 



|F £ (r,£)|<C, 



< C^ 2 " 1 



we verify - discussing the cases sgn £i = sgn £2 and sgn £1 / sgn £2 separately - the estimate 

\F £ (t, e + C2) - F e (r, £ + Ci) I < Ce 72 " 1 K2 - Cil < C K2 - Ci| 72 (3-9) 

for all r G /, £ G eZ, and all £i, £2 £ [ — e, +e]- I n particular, setting £1 = and taking the supremum 
over £ and £2 we obtain the convergence estimate 

\\F £ -F £ \\ L oo {n) <Ce^. 

We next show that F £ is Holder continuous with respect to £. To this end, let r G / and £i,£2 G M 
be given and denote by £« the natural projection of £j to eZ, that means |£i — £«| < e/2. For 
|£i — £2! < e/2 we have £1, £2 G [£1— e, £i+e] and are hence done by (3.9). In the case of |£2 — £i| > e /2 
we combine our assumptions of F £ with (3.9) and the triangle inequality to obtain 

\f £ (t, 6) - f £ (t, ei)| < c |£ 2 - £i| 72 + Ce^ 2 -^ |6 - 6| 72 + |6 - £ 2 | 72 ) 

<C|6-6| 72 + Ce^ 2 
and the desired estimate follows with 

|a - 6| < 16 - ai + 2 • e/2 < 3 16 - ai . 

The claim of the Lemma is now a consequence of the convergence estimate, the spatial and temporal 
Holder estimates for F £ , and the Arzela-Ascoli theorem applied to F £ . □ 

Since the functions Q £ are obtained by solving the discrete heat equation with macroscopic initial 
data, they converge as e — > to a smooth solution of the macroscopic heat equation d T Q = d?Q. 
This is not surprising and can be proven in many different ways. For our purposes, it is sufficient 
to observe that strong compactness is provided by combining Lemma 3.10 with the following Holder 
estimates. 
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Lemma 3.11 (Holder estimates for Q £ ). There exists a constant C independent of e such that 

\Qe(T2, 6) - Qe(n, £i)| < C{ \t 2 - nl 1 / 2 + |& - &| ) 
holds with £1,62 G ^Z aric ^ < Ti, T2 < Tfi n /or a// < e < £*. 
Proof. From Corollary 3.7 and |<7o(£)l — C£~ we derive 

*t 2 



\ gj (t 2 ) - ^(ti)| < / 2 |^(t)| dt < ae 2 (t 2 -h) + C/3e( v / t 2 "- v 7 ^), 



as well as 



J2-1 
|?»(*i) - %i(ti)| < X] l V +9i( i i)| < ae(j 2 -ji) 

3=31 

for all ji, j'2 G Z with ji < j'2 and all < t\ < £2 < 00. Setting jj = £j/e and ti = Tj/e 2 we therefore 
get 



\Qs(T2, &) - Qe(n, Ci)| < <*(& - &) + a(r 2 - n) + 0(v^2 - v 7 ^), 
and the claim follows since < n < T2 < Tfi n implies |t2 — ri| + (v/^2 — v/ri) < C\Jt 2 — t\. D 



A crucial part of our analysis is to establish strong L°°-compactness of the functions R £ . The 
main difficulty is that R £ equals the sum of K £ shifted and delayed versions of the discrete heat kernel 
producing a temporal discontinuity at any of the phase transition times rj" £ , . . . , r^- . In order to 
control the impact of all these jumps in time we split 

R £ (T, = Rl, e(T, O + Ri, e(T,0 
with 



R h£ (r, := -2^fl- e (r-r fc * e , C-efc) 



fc=i 



and i?2, e := Re — Ri,e- Here the function H £ : R 2 

H £ (t, := < 



( for r < 0, 

7" 

-p- G £ (d*e, for < r < d*£, 

k G £ (r, for r > d*e 



can be regarded as a temporally regularized version of G £ , where the latter represents the discrete heat 
kernel in macroscopic coordinates according to G £ (e 2 t, ej) = gj(t). In particular, H £ is continuous 
with respect to r G M and differs from G £ for < r < d*£ only. 

The function R 2 , e contains all the temporal jumps caused by the phase transitions and can 
therefore not be compact L°°(f2). The key observation, however, is that R 2 ,e is still uniformly 
bounded and converges as e — > to in L s (il) for any 1 < s < 00. The macroscopic limit of 
(R £ ) 0<£<£ is therefore completely determined by the family (i?i )£ )„ < , where each function i?i i£ 
is continuous with respect to r and satisfies 

Q £ {r* k £ , ek) + i?i, £ (r fe * £ , ek) = lim P £ (r, ek) = 1 (3.10) 

T /~ T k, e 

for all k > 1 with t^ £ < rfl n thanks to limj ;*j* Pfc(t) = 1- 
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Lemma 3.12 (bounds for B,2,e)- There exists a constant C independent of e such that 

II #2, e||L°° (ft) < C, H^ellL^fJ) < Ce 

holds for all < e < e* . 

Proof. By construction and Corollary 3.5 we have 

SUppi?2,e C (ll, e U ■■■\JlK e ,e) X K, Ifc, £ := [Tfc j£ , t£ £ + d*e], 

and the intervals ifc e are pairwise disjoint thanks to (3.8). In particular, using 

R2,e{T~k £ + o~, je) = 2H £ (a, ej — ek) — 2G £ (a, ej — ek) for all j € Z, <r < d*£, 
we estimate 

|^2, £ «e + CT > J e )\ < 2 (G £ (0, ej - ek) + G e (d,e, ej - efc)) < 2, 
see Lemma A. 2, as well as 

/ \R2,e{r* k , e + <r, £)| d£ < 2 £ ^ (<&-k(0) + 5j -_ fc (d,/e)) = 4e 

thanks to YsjGz9j(t) = 1 for all t. The first estimate implies ||i?2,£||L°°(n) < 2, whereas the second 
gives rise to 

^i rd t e r 

\\R2,e\\ L i(n) < Y\ / / |^ 2 , e (r fc % + a, £)| d£dcr<2if £ d*£ 2 <T fln £, 

fc=1 JO JR 

where we used (3.8) again. D 

It remains to establish L°°-compactness results for R\ y£ . To this end we next derive a further 
auxiliary result concerning the Holder continuity of H £ . 

Lemma 3.13 (Holder estimates for H £ ). For each < 7 < 1 there exists a constant C independent 
of e such that 

I tt 1 c \ tj 1 t \\ s- n I \ T 2- Ti\ S2 -?1 \ 

\H £ {t2, 6) - H e (n, £1) < Ce[ — — — - + 



max{d ! „e,Ti}')'+ 1 / 2 maxjd^Ti} 3 / 4 / 

holds with £1, £2 £ £^ and < n < T2 < Tfi n /or a// < e < £*. 

Proof. Suppose at first that d*e < t\ < t 2 . Thanks to the temporal Holder estimates for the discrete 
heat kernel, see Lemma A. 3, we find 

\h £ (t 2 , & - H £ (n, b)\ < V g J/ 2 =^ 2 \T2-nr 

\e 2 ) T± 

for some constant C independent of e and £2- Similarly, Lemma A. 3 also ensures that 

1/2 



C 



e e 



Ce 



\H e ( Tl , 6) - H £ (n, a)| < ;_ 3/4 - -^ 16 - £i| 1/2 



©f 



r- 
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Now suppose that < n < r 2 < d*e. We then estimate 

\h £ (t 2 , 6) - H e { n , 6)| < g6/£ . (d ; /£) (r 2 - n) < - §7 ^— (r 2 - n) 

a*£ ^/ £ i/2 

C Ce 

< -WJ7, ~{d*e) ~ 7 |t 2 - Tl| 7 < — - -^ |T2 - n| 7 

d^e 1 / 2 (d,e) 7+1/2 



as well as 



C 

i Tl i i 

\H e {n, 6) - H £ (n, 6)| = -j—\9&/e{d*/e) -9^/e{d*/e)\ < 



e e 



1/2 

c £ |6-6l 1/2 



4 £ '^ £ — ^i/«V-/^,-: ^3/4 (4£) 3/4 

£ 

In summary, we have established the desired estimates in the special cases < n < r 2 < d*£ 
or d*£ < ti < r 2 . All other cases can be easily be traced back to these cases using the triangle 
inequality. □ 

We are now able to prove our main technical result in this section. 

Lemma 3.14 (Holder estimates for R\ )£ ). F° r each < 7 < 1/2 there exists a constant C indepen- 
dent of e such that 



\RiAt2, 6) - RiAn, 6)| < C( \t 2 - rip + |6 - 6 

holds with 6)6 £ £ Z an d < ri < r 2 < r-a n for all < £ < e* . 



1/2" 



Proof. It is sufficient to proof the assertions in time and space separately. 

Holder continuity with respect to £: Let < r < rfi n and 6^6 S M be given. Then there exists 
m £ £ Z such that 

'm s ,e ^ ' — 'm e +l,£ 

and (3.8) ensure that em e < r/(2d*) as well as 

r — rjt £ > 2d*e{m £ — k) for all fe = 1, . . . , m e . 
In particular, we have 

^( r - T M'0 =0 for a11 fc>m ^ ? GM ' 
so Lemma 3.13 yields 

\Ri,e(r, 6) - Ri,e(r, 6)| < 2 J] |^(t - < e > 6) - H £ (t - r fe %, 6)| < CS 1 , |6 - £i| 1/2 

fe=i 

with 

„ v^ e <- Y^ < I < r 



where we used the Riemann sum approximation of the integral as well as the monotonicity of the 
integrand. 

Holder continuity with respect to r: Now let £ £ M and < ri < r 2 < T^ n be fixed, and choose 
m £ , n £ £ Z such that 

T m s ,e <T 1- T m e + l,e> T n £ ,e ^ T 2 - r n e + l,e- 
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This gives 

em £ < n/(2d*), e{n £ - m £ ) < (t 2 - n)/(2d*) 

as well as 

\Ri,e(r 2 , f) - i?i, e (ri, 0| < 2A £ + 2Y £ , 
where 

^■■=Y.\ H ^-rl £ ,i)-H £ {r 1 -Tl £ ,Ol Y £ := £ |F £ (r 2 - r fc %, £) |. 

fe=l fc=m e +l 

Similar to the above, we deduce that 

v ^^ Ce|r 2 -ri| 7 ^ / 1 1 p da V .. ^ . |7 

*< ^ E maxR£ , T _ r * e}7+V2 ^U + ^ io ^J^-^^-^i 7 ' 

whereas Y^ can be estimated by 

A;=m £ + 1 fc=m e +l 

2d*e J T1 2d*e J Q 

z C [ 2T1 % = c ^~ ri ' 1/2 <C\T2-n\\ 

see Lemma A. 2 and (3.8). □ 

We conclude this section by showing that both the L 1 -norm of i?i ;£ (r, •) and the L 2 -norm of 
V e i?i e (r, •) are bounded uniformly with respect to r and e, where 

VeF(f) ;= F(f + ,)- fK ) 

is the discrete spatial gradient of a function F defined on R. 

Lemma 3.15 (Lebesgue bounds for R\ £ and its discrete gradient). There exists a constant C such 
that 

ll-frl.ellLo^J;!. 1 ^)) + II V £ i?l, £ ||L°°(/ ; |_2(R)) < C 

holds for all < e < e* . 

Proof. For r > ed*, the properties of the discrete heat kernel, see Lemma A. 2, imply 

f H £ (r, Ode = eJXr/e 2 ) = e, / (v £ tf £ (r, £)) V = e" 1 £ (v+^(r/e 2 )) 2 < ^, 

and combining this with the definition of H £ for all r G R we find 

{0 for r < 0, 

£ V4 for < r < d,e, 

er -3/4 f Qr T > ^ £ 

From the first estimate we infer that 

\\Rl,e(T, -)IIlI(K) - £Ks -^f 
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holds for all r € I. We next fix t £ I and choose an integer m £ such that r^ e < r < r^ +1 e . As 
in the first part of the proof of Lemma 3.14, we estimate 



|V £J Ri, e (r, -)lb(R) <£||V £J ff ( 



4 T_r fc,£> -JllL2( 



fc=l 

771^ — 1 



* C£V4 + E ,„,,., ,,,3/4 = ^ 1/4 + 1 



tZi (2kd*s(m e - k)Y 
and the proof is complete. □ 



3.4 Convergence results and verification of limit dynamics 

In view of the compactness results in Lemmas 3.9, 3.11, and 3.14, we may select a subsequence of 
e — > 0, which we do not relabel, such that 

£ -> T in L°°(I), Q £ -). Q in L°°(fi), P M -> i? in L°°(fi). (3.11) 

As i?2,e —5- in L s (f2) for any 1 < s < oo by Lemma 3.12, we find 

P £ = Q £ + R h£ + R 2 , £ ^Q + R = P inLf oc (tt), (3.12) 

that means, P is the limit of P £ in Lj' oc (il) and the limit of Q £ + R\^ £ in L°°(Q). Moreover, convergence 
of (Q £ ) £ implies convergence of the initial data 

P £ (0,-) = Q £ (0,-)^Q(0,-) = P(0,-) inL°°(R). (3.13) 

Theorem 3.16 (limit dynamics along subsequences). Any limit (P,Q,R,(,*) satisfies: 

1. S* := {(t,£*(t)) : t £ 1} is a Lipschitz curve in $7; the functions Q, R, and P = Q + R 
are bounded and locally Holder continuous in £1; furthermore, R € L°°(/;L 1 (M)) and d^R 6 
L°°(/;L 2 (M)). 

2. Q is a solution of the heat equation in Q with initial data P(0, •). 

3. (P, £*) is a distributional solution of 

d T P = df P in ft \ S*, \[P}\ = and 2±C = |[%P| on E* (3.14) 

with initial data (P(0, •),£*(())) attained in L°°(M) x IR. Moreover, we have 

P(r,0>-1 for all{T,i)£tt, 

p(r,o<i */e>r(r), 

which implies P G [—1,1] oraH*, and t/ie movement of the interface is determined by 

£fC(. T ) > /or a/mosi a// t£I, 

Ae*(r) = i/P(r,r(r))^l- 

Remark 3.17. Being a distributional solution of (3.14) means 

-/ (P + fi)d T i;d^dT= / PdfV>d£dr, ^C»((0,r fil ,)xE), (3.15) 

JO JR JO JR 

where /i(r, £) = sgn(£*(r) — £). In the following we will use £* and \i interchangeably to represent 
the solution, whichever is more convenient. 
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Proof of Theorem 3.16. The continuity properties of £* and Q are immediate consequences of 
Lemma 3.9 and Lemma 3.11; Lemma 3.9 also yields g^£* > 0. Holder continuity of P and the 
claims for R follow from P = Q + R and the bounds on R\ )£ proved in Lemmas 3.14 and 3.15. 



Setting h e {t,0 = sgn(£*(r) 
distributional form as 



£) = sgn U £ (r, £) , we write the equation for U £ = P £ + fi £ 



in 



Tfin 



(P e +/i e )^dedr 



Tfin 



pA^d^dr, yeq°((o,T fin ) 



and deduce (3.15) in the limit £ — )■ 0. Similarly, Q solves the heat equation, and both P and Q attain 
their initial data in L°°(]R) due to (3.13) and continuity of Q. 

By construction, the discrete solutions satisfy P £ > — 1 in Q and P £ < 1 in {(r, £) G fi : £ > £* (r)} 
for all £ > 0, and in the limit e — > we obtain the corresponding inequalities for P and £*. In 
particular, we have P(r, £*(t)) £ [ — 1> 1] f° r an t G I. 

It remains to check that P(r, £*(r)) < 1 implies ^£*(t) = 0. To this end, let (t,£) G H* and 
(5 > be given such that P(f,£) = 1 — 25, and suppose at first that £ < £*(rfi n ). For any e choose 
£ e G eL such that £ e < £ < £ e + £ and denote by f £ the phase transition time corresponding to £ e . 
See the left panel of Figure 3.4 for an illustration and notice that f £ < Tfi n because otherwise the 
interface position would be maximal via £ = £*(Tg n ). Holder continuity of P now implies 

P(T £ ,& -1 + 25 < \P(f £ , &) - P(f , 6| < C(|f £ - f P + £ t) 



for some exponent < 7 < 1, while uniform convergence of Q £ + Pi, 
Ri,e)(f e ,£ e ) = l, see (3.10), yield 

P(f e ,&)>l-o(l) e _> - 



P as £ — )■ and (Q e + 



We thus find 5 < C\t £ — r| 7 for all sufficiently small £ > and may select a subsequence of e — > 
such that r e — > tq and fo < f . The uniform convergence £* — > £* implies 



r(r ) = hW*(r £ 







lim £ F 

e^0 



e(f) 



and gj£* > ensures that £* is constant in [to,t]. By a similar argument using £ e = £ e + £ and 
the corresponding phase transition time f £ we finally conclude that f is a regular point of £* and 
g^£*(r) = 0. Moreover, in case of f < rg n and £ = £*(rfi n ) we find that £* is constant on [f, Tfin], 
and for f = Tfi n we can repeat the above reasoning after enlarging the time interval slightly beyond 

Tfin- □ 
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Figure 3.4: Left. Illustration of the key argument in the proof of Theorem 3.16: By construction, we have 
£ = £*(t) = £*(to) and -§^C > 0, so £* is constant on [fo,f]. Right. Approximation s n of the sign function 
used in the proof of Theorem 3.18 for the case £J > Q\ f° r £1 ^ £2 one has to redefine s n such that s n (0) = — 1. 

We complement Theorem 3.16 with a uniqueness result by adapting some techniques for hysteresis 
problems from [Hil89, Vis06]. 

Theorem 3.18 (well-posedness of the limit problem). The solution to the limit problem in Theo- 
rem 3.16 is uniquely determined by the initial data P(0, -)and £*(0). 
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Proof. Given two solutions (P±,Qi, Ri,Q) and (P 2 , Q2,R2,^) wrfcn initial data Pi(0, •) = P2(0, •) 
and £i(0) = £ 2 (0), we se t P = P\ — P2 and Jl = ^1 — ^2, where /^(t, £) = sgn(£*(r) — £) as in Remark 
3.17. In order to show P = p, = we follow the strategy of [Hil89, Theorem 5], that means we first 
establish sufficient regularity in time and derive afterwards an L 1 -contraction inequality by testing 
the equation for P with sgn P. 

Regularity in time: Standard uniqueness results for the heat equation imply Q\ = Q2, and we 
find P = Ri -R 2 G L°°(I; L^R)) and <%P G L°°(/; L 2 (M)) in addition to boundedness and continuity 
Furthermore, J1(t, £) is bounded, and it is nonzero only if £ lies between £J(r) and £f(r). We thus 
conclude 

PG L 00 (/;L 1 (R))nL 00 (/;H 1 (R)) and js G L°°(I; L 2 (M)). (3.16) 

In view of (3.16), we may integrate by parts after subtracting the equations for (Pi, m) and (P2, // 2 ) 
from each other, which gives 

/ n ^(r) / (P + AZ)(r,Or?(Od£dr= / ^(r) / %P(r, £) %(£) d£ dr 

jo Jm Jo Jr 

for all (p G C£°(P), 77 G C£°(R), and by density also for all (f> G Hj(I), n G H^M). Consequently, 

4- / {P(r, + A(r, £)) r/(£) d£ = - / %P(r, £) %(£) d£ (3.17) 

a r ii Jr 

for all t £ I. A direct computation shows that 

-f / Hr,0v(0^ = MCi(r)) ££(r) - 2r ? (£ 2 *(r)) ££(t) (3.18) 

ar JR 

for all r G / where g^^(r) and ^£ 2 (r) are defined, and the right hand side of (3.18) can easily 
be bounded by C 1 1 77 1 1 h 1 (]R) 5 where the constant C depends on ||£?||w 1 .°°(/)) J = 1)2. Thus, <9 T P(r, •) 
exists in H 1 (M.) and 



[d T P(T,-), V ) 



£ [ P(t,0v(04$ < C (||^^|| L oo (J;L2(M)) + 1) ||/7||h1( 

Jm. 



By standard embedding results, see for instance [Eva98, Thm. 3 in Sec. 5.9] and note that / is 
bounded, d T P G L°°(J; H" 1 ^)) and P G L°°(J; H^K)) imply P G C(I;L 2 (R)), and with (3.16) we 
conclude PG C(I;L 1 (R)). 

Contraction inequality: Given r G [tl,T2], where < n < r 2 < Tfi n such that Q > Q m [ti,t 2 ], 
we approximate the sign function by 



s n (z) 



max(— 1, min(l, 1 + nz)) if \z\ < n, 

n + sgnz— \z\ if n < \z\ < n + 1, 

otherwise, 



see Figure 3.4 for an illustration. In what follows we suppose n > ||P||oo and consider n = 
s n (P(r, •)) G H^R) in (3.17)-(3.18). Due to P(t,Ci(t)) > if ^*(t) > and s„(0) = 1, we 
then find 

^i>))^*( r ) > lCi(r), 

while s n < 1 and ^£ 2 (t) — imply 

-m(T))±e2(r) > -;fc£(T). 

Hence, we obtain the Hilpert estimate 

»K£(t))&£(t) - ri{^ T ))^{ T ) > &£(r) - &&(?), 
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and as moreover the right hand side of (3.17) is nonpositive, we infer 

(d T P(r, •), s n (P(T, •))> + 2 (££(r) - £C 2 (r)) < 0. (3.19) 

Using 

(d T P(T,-),s n (P(r,-))) = A / S n (P(r,0)d& 

where <S^(z) = s n (z) and 5^(0) = 0, we next integrate (3.19) from n to T2 and arrive at 

[ S n (P(T 2 ,0)d£ + 2\e i (T 2 )-e2(T 2 )\< /5 n (P(ri,0)de + 2|ej(ri)-e|(ri)|. 

By construction, S n (P(r, •)) converges to |P(r, -)| in L 1 (M) as n — )■ oo, and passing to the limit yields 
the desired inequality 

||P(r 2 ,.)|| L1 (R) + 2|^(r2)-^(r 2 )| < ||P(n, .)|| L iCR) + 2 |£(n) - 3(n)| (3.20) 

in the case of £J > £ 2 m [tii - ^]- Moreover, for £* < Q i n [ T i) T 2] we derive (3.20) by repeating the 
above arguments with s n (z) = max(— l,min(l, — 1 + nz)) for \z\ < n, which satisfies s n (0) = — 1. 
Combining both cases and continuity of ||P(t, •)IIl 1 (r) we finally obtain 

ii^iIl«.(j;li(r)) + 2|i£ - ? 2 *iiloo ( 7) < iip(o, -)\\ L i m + 2 ier(o) - &m , 

so uniqueness follows from P(0, •) = and £f(0) = £1(0)- D 

As a consequence of Theorems 3.16 and 3.18 we obtain the following approximation result. 

Corollary 3.19 (uniqueness and improved convergence). If P £ (0, •) — > P(0, •) in L°°(IR) as e — > 0, 
the limit (P,Q,R,£*) in Theorem 3.16 is unique and the convergence (3.11)-(3.12) holds along the 
whole family e — >■ 0. 

A The discrete heat kernel 

Denoting by ~ the Fourier transform with respect to the discrete spatial variable j, that is 

g(t,k) = Y,9j(t)e-* kj , k€[-n,n), 

the initial value problem (3.4) for the discrete heat kernel transforms into 

d t g(t, k) = p(k)g(t, k), £(0, k) = 1, (A.l) 

where p is the Fourier symbol of the negative discrete Laplacian, i. e. 

p(k) = 2 - e~ lk - e +lk = 2(1 - cosfc). 
Solving the parametrized ODE (A.l) and applying the inverse Fourier transform we find 

9j(t) = ^J ^(t, k)e lkj dk=±-J \xp (-p(fc)t) cos (jk) 6k. (A.2) 

Lemma A.l (monotonicity and convexity properties for j = 0). 

1. 1 1—?- go(t) is strictly positive, strictly decreasing, and strictly convex. 

2. 1 t— > L go(s) ds is strictly positive, strictly increasing, and strictly concave. 
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3. 1 1— >■ go (t) is strictly negative, strictly increasing, and strictly concave. 
Proof. The representation formula (A. 2) implies 

go(t) > 0, g {t) < 0, go(t) > 0, g (t) < 

for all t > 0, so all assertions follow immediately. □ 

Employing standard methods from asymptotic analysis one finds 

as well as asymptotic laws for the long-time behavior of any discrete moment. For our considerations 
in §3, however, the following rather rough estimates are sufficient. 

Lemma A. 2 (temporal decay properties). There exist positive constants c and C such that 

<<?,(*) <3oM<C(l + tr 1/2 , (A.3) 

| A 9j (t)\ = \ 9j (t)\ < -g (t) < C(l + 1)- 3/2 (A.4) 

and 

gv(t)>c{i + ty l l 2 (A.5) 

hold for all j G Z and all t > 0. Moreover, we have 

J>(t) = 1, £ |V +5 ,(i)| 2 < C(l + tr 3 / 2 (A.6) 

/or all t > and some constant C . 
Proof. For £ > 1 we observe that 

-k 2 < ^k 2 < p{k) < k 2 for all fc€[-7r,7r], 



1 /"+7T i r+oo 

\9j(t)\ < go(t) < — / exp ( - fc 2 t/4) dfc < — = / exp ( - fc 2 ) dfc 



7T 



and this implies 

»+oo 

PY n r _ k 2 \ Ah = 

as well as 

„ W > JL £* exp ( - Ift) d fc = ^ £*^ exp ( - e) d fc > -L^ £' exp ( - *») d*. 

Combining these estimates with < <?o(l) < 5o(i) < 5o(0) = 1 for < £ < 1 we readily obtain (A.3) 
and (A.4). Moreover, for t > 1 we estimate 



A 5j(*)l = \<Jj 



[t)\ < ^ y + " P(fc) exp ( - p(Ar)t) dfe = -5o(t) < ^ y + " A: 2 exp ( - A: 2 t/4) dk < ^ 



and this provides (A.5) due to < —go(t) < — <?o(0) = — A^o(0) = 2 for all t. The discrete heat 
equation (3.4) further ensures conservation of mass via Yljez9j(t) = J2jez9j(®) = -*-• ^- n particular, 
using discrete integration by parts as well as Holder's inequality for series we find 

^(V+^(t)) 2 = -J2 gj (t)A 9j (t) < IIA^Hoo 

jez iez 

which implies the estimate in (A.6) thanks to (A.4). □ 
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A further key ingredient to our convergence proof in §3 are the following time-dependent Holder 
estimates for Qj{t). 

Lemma A. 3 (longtime behavior of spatial and temporal Holder constants). For each 7 £ [0, 1] there 
exists a constant C 7 such that 



M*2)~gj(*l)| < ,, f| : , ,-,-J 2 

t 2 >*i 1*2 — h\ 



sup j- py L < C 7 (l + £ 



/joWs /or a// j £ 7* and all t\ > 0. Moreover, there exists a constant C such that 

|gj 2 (*)-gji(*)| . ^ n , .v-3/4 



sup '^ w ^^ <C(l + t)- 
1 



J1J26Z |i2-i jl " 




/loWs /or all t > 0. 

Proof. Let j £ Z and < £1 < £2 be fixed. Thanks to (A. 4) we estimate 

\ 9j (t 2 ) - gj (h)\ < 1 2 \ gj (t)\ dt < c f 2 (1 + 1)- 3 / 2 dt = c|(i + ti)- i/2 - (1 + t 2 r 1/2 |, 

and writing 1 + £2 = s(l + £1) with s > 1 we get 

\ 9l (t2) - g 3 (ti)\ C/ 7 (s) _ 1 g V 2 - 1 

|£ 2 -£i| 7 - (1 + ^)^/2' M«;- s i/2 (s _ ir - 

We readily check that the function / 7 is bounded on [1, 00), so the first claim follows by taking the 
supremum over s > 1. 

Now let £ > and Ji, J2 & %• with j 2 > f\ be arbitrary. By Holder's inequality for series we then 
find 

A2-1 \ 2 fh-i 

W*)-fti(*)| 2 = El v +»(*)l ^ £ 

V=ii / \j=h 

< \T,\ V +9j(t)\ 2 )(j2-Ji), 

\J6Z / 

and the second assertion follows from (A. 6). D 
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